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Abstract 



Long range Coulomb interaction between the edges of a Hall bar changes 
the nature of the gapless edge excitations. Instead of independent modes prop- 
agating in opposite directions on each edge as expected for a short range inter- 
action one finds elementary excitations living simultaneously on both edges, 
i.e. composed of correlated density waves propagating in the same direction 
on opposite edges. We discuss the microscopic features of this Coulomb drag 
of excitations in the fractional quantum Hall regime within the framework of 
the bosonic Chern-Simons Landau-Ginzburg theory. The dispersion law of 
these novel excitations is non linear and depends on the distance between the 
edges as well as on the current that flows through the sample. The latter de- 
pendence indicates a possibility of parametric excitation of these modes. The 
bulk distributions of the density and currents of the edge excitations differ 

significantly for short and long range interactions. 
PACS numbers: 73.40.Hm, 72.15.Nj 
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The Integer and Fractional Quantum Hall Effects occur in a two dimensional electron 
gas placed in a strong perpendicular magnetic field. Interesting phenomena can take place 
at the edges of this gas. The simplest is related to the edge currents which are instrumental 
in our understanding of the transport properties of the QHE ]1|-|J. The character of the 
excitations in the Quantum Hall system is also effected by the presence of the edges. While 
the bulk excitations have a finite gap the excitations which are localized near the edges are 
found to be gapless. They are chiral, that is moving only in one direction along the edge, and 
can be described by the chiral Luttinger model f| provided one disregards their coupling to 
other modes. Recently the possibility to probe the structure of the fractional Hall states by 
studying the tunneling between edges attracted much attention ||[7]]. 

In this paper we study the edge excitations within the framework of a bosonic Chern- 
Simons Landau- Ginsburg (CSLG) theory. This theory || was proposed as a useful sup- 
plement and extension of Laughlin's fully microscopic many body theory of the FQHE. 
Recently, the edge effects were studied using this theory || and Wen's results were red- 
erived. Our goal is to investigate the effects of the inter-edge interactions. In pi novel 



effects were predicted due to this interaction which are related to the redefinition of the gap- 
less modes into a Bogolubov transformed combination of modes living on both edges and 
propagating in the same direction. Our analysis goes one step further and represents the 
microscopic picture of this phenomenon. We use the realistic Coulomb interaction between 
the edges and derive the modified modes and their density and current distributions. These 
clearly indicate that one deals with the Coulomb drag of excitations on one edge by the 
other. We derive a non linear dispersion relation for the modified excitation modes and find 
that it depends on the width of the Hall bar and on the total current which flows through 
it. We point out that the latter dependence opens a possibility of a parametric excitation 
of the novel modes. 

The starting point of our analysis is the mean field equations derived by minimizing 
the CSLG action for the fractional quantum Hall effect ||. It will be convenient to write 
these equations in a hydrodynamic form using the polar decomposition of the bosonic field 



<p = y/pe -19 and introducing the velocity fields Vi = ^Qfi + ^( a i + A*) with the vector 
potential A taken such as to produce a constant magnetic field along the z direction, 



u<%% = 1 - P ( la ) 



- 9 t t)j - d^o = eijpvj (lb) 
dtp = -di(pvi) (lc) 



\[f> 

Here and throughout the paper length is measured in units of the magnetic length I 



d t 9 = -\v\ + + a - V(x) - p -JLJu(t- r')p(r')dV . (Id) 



fic/eB, time in units of inverse cyclotron frequency u c = eB/mc, energy is normalized by 
hu c and the density by its bulk value p = uB/ipo where v = 1/ (2n + 1) is the filling factor. 
All the quantities apart of 9 and ao are gauge invariant and the latter appear in the gauge 
invariant combination d t 9 — a . In order to correctly reproduce the energetics of the nonin- 



teracting limit we have followed Ref. |TTJ and included a 5-function type repulsive force (the 
term before the last in ([Td|) ) with a strength of 2nTi 2 / vm. Our confining potential is assumed 
to rise fast enough to avoid the occurrence of an alternating sequence of compressible and 
incompressible strips along the edge 0,|J. Apart from this requirement the detailed nature 
of V(x) is of no importance for our applications. 

Our aim is to obtain the edge excitations as the RPA modes of the theory , i.e. the 
eigenmodes of the above equations linearized around a static solution. We first consider 
a single edge. The translation invariance in the y direction suggests looking for a static 
solution in the form: 

9 = -x y - p(x )t , p = p(x) , v x = , v y = v y (x) , a = a Q (x) . (2) 

We assume for definiteness an infinitely high wall situated at x < and accordingly set the 
density to zero at the wall. The only gauge freedom which preserves the form of the above 
solution is adding arbitrary constants to x$ and p. We fix this freedom by requiring that the 
statistical potentials a , vanish at x = and choosing A = (0, Bx, 0). This assures that 
solutions with different values of xq and p are not gauge transforms of each other. Formally 



xq is the conserved momentum along the edge and its physical effect on the condensate <p is 
similar to the guiding center coordinate of Landau levels - changing its value translates </> in 
the x-direction. For a fixed xq the value of p is determined by requiring that p approaches 
p far from the edge. It can be shown that p is the energy which is needed in order to add a 
particle to the edge. Representative examples of the density and current density profiles for 
solutions @) in the case of short range interactions U = (2irh 2 \ s /ism)5(r — r') are shown in 
Fig. |l|. One finds a one parameter family of static solutions depending on xq and differing 
by the density of particles at the edge. With our boundary conditions the velocity at the 
edge is v y (0) = —xq- Using this and the fact that the velocity falls to zero in the bulk we 
integrate Eq. flla|) and find that — x is the excess charge per unit length along the edge 
relative to a step-like constant density profile. 

For a Coulomb interaction U = 2-Kv^ 1 \ c huo c / | r — r' | with A c = ue 2 /2irel Hu c and a 
constant neutralizing background extending up to the wall one obtains distributions similar 
to Fig. |I[ The essential difference between the two cases is found in the large x behavior 
of the solutions. By a simple analysis of Eqs. ( JIa| - |Tcl| ) one finds that for the short range 
interaction all the relevant quantities decay exponentially to their bulk values with a decay 
constant a = ^2(X S + 1) — [\ s + l) 2 — 1. For the Coulomb interaction the density and 
the velocity behave at large distances according to 1 — 2\ c xo/x 2 and — 2\ c xo/x respectively. 

Even more pronounced difference is found in the dependence of p on xq which is displayed 
in Fig. ^. For the short range interaction the chemical potential tends to |(A S + 1) for 
large positive values of xq and increases as the charges are pushed against the wall, i.e. 
for decreasing xq. In the case of the Coulomb interaction, assuming that the length of the 
sample L is very large compared to xq, one finds that the leading contribution to p is given by 
— 2X c xo In L. This term corresponds to the electrostatic charging energy of the excess charge 
(—Xq) had it been uniformly distributed over a strip of width I along the edge. Its change of 
sign at Xo = reflects the tendency of the system to remain neutral. The remaining part of 
the chemical potential p = p + 2A c Xo In L is due to the kinetic energy and the deviation of 
the excess charge distribution from that of a uniform strip. As one can observe from Fig. |2| 
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the dependence of p on xq varies considerably as the strength of the interaction is changed. 
Linearizing Eqs. (p^aj-p^ around one of the solutions (|2|) we obtain the RPA equations 



€ijdi5vj = —5p (3a) 
d t di59 - d t 5vi - di5a = e iy v y 5p + e^pdvj (3b) 



d t Sp = -v y dy5p - di(p5vi) (3c) 

djSp 
Ap 



d t 56 = -v y 5v y + f(p, 5p) + ^ + 5a -5p-^ f U(r - r')5p(r')d 2 r' . (3d) 

4*0 271 J 



where / denotes the part of the linearized "quantum pressure" term, i.e. the second term 



on the r.h.s. of (|TdT ), containing x derivatives. 



Our main method of solving Eqs. (l3a| - |3d| ) will be to use the continuous one parameter 
family of static solutions described above in order to find the gapless branch of the RPA 
eigenmodes. The derivative of the solution (0) with respect to Xq constitute a static solution 
of the RPA equations (g|-0). We denote it by 56 {s) , 5p (s) (a;) , 5v y s \x) , 5a { s \x). Motivated 
by this observation and concentrating on the long wavelength limit we set 

5p = 5p^ s \x)cos[ky — oot) 

5v y = 5v y s \x)cos(ky - out) (4) 
5a = Sclq (x)cos(ky — oot) 
56 = ——sin(ky — ut) . 

Clearly these density and velocity distributions are concentrated along the edge. With the 
expression for u given below also the gauge invariant combination d t S6 — 5a vanishes 
far away from the edge. Inserting these functions into Eq. ( pd|) and neglecting the term 
proportional to k 2 we find that they indeed represent a solution of the equation provided we 
choose properly the dispersion relation u = u(k). In the case of the short range interaction 
we find directly 

" = -£H*< 1+ ">- (5) 
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where we have restored the units of dimensions. In this expression we have included the 
second order term in k which can be obtained after multiplying the first three fields in 
(|j) by 1 + kl. In the case of the Coulomb interaction after the substitution of 5p the last 
term in Eq. fl3cf ) becomes —2X c cos(ky — cut) / °° dx'8p( s \x')KQ(\k\\x — x'\) where Kq(x) is the 
modified Bessel function. Working in the region x <C k~ x and assuming that k~ l is much 
larger then the width of the region where Sp^ is appreciable we can use the approximation 
i^od^H^ — ~ ln(2e~ 7 /|/;;| \x — x'\) where 7 is the Euler constant. As a result we find 
that Eq. (fkj) is satisfied to first order in k when the dispersion law is modified by an extra 



logarithmic term, cf., Refs |T^JT3|] 



Note that p rather then p (cf., Fig 0) enters this expression. 

Integrating the linearized continuity equation fl3c]) we find that the above solution should 
be supplemented by a velocity field in the x direction: 

5v x = —{ijj5v^ — k5a^ + uj)sin(ky — oot) (7) 
p 

which corresponds to a current pSv x vanishing on the wall as well as far away from it. One 
can also check that the other two equations in the linearized set are satisfied up to first order 
in k. 

For the short range interaction the linearity of the dispersion relation and the fact that 
dp/dxo < imply the chirality of the waves. The wave velocity is determined by the 
compressibility of the edge k = (px^ dp/dxo)" 1 which is controlled by the wall potential. 
It vanishes when the edge of the static solution is far away from the wall (i.e. at large 
positive xq) for which case u ~ k 3 in agreement with Ref. fL2|| . In the Coulomb case the 



long wavelength limit of the dispersion curve is dominated by the logarithmic part so that 
although dp/dxo can have either sign the excitations are chiral as before. For distances 
which are larger than k~ x the Bessel function Kq decays exponentially and we expect a 
crossover in the excitation profiles from the Coulomb power law behavior to an exponential 
decay. 



We go on now to consider the case of a wide and long Hall bar defined by hard walls at 
x = ±-y. According to our strategy we first show that again a family of static solutions of 
Eqs. ( p^aj - p^D exists. This will now be a two parameter family which we find by appropriately 
gluing together two solutions of the single wall case. Assuming that the two solutions 
correspond to the values of xq and \x which are (x\, fix) and (22, ^2) respectively and denoting 
the solutions by superscripts 1 and 2 we set 

9 = —x\y — nit (short range) , 9 = —X\y — (/ii — 2\ C X2 In (Coulomb) 
p = p (D( x + y)8(-x) + p (2) (y - x)e(x) 

v x = (8) 

v y = v^{x + —)-v^{^--x) 

a (x) = a^\x + y) + 4 2) (^ - x) - 0%\W) 

where G is the step function and where we assumed W <ti L. The gauge fixing in the present 
case is achieved through the requirements that vanishes on the left wall and a y equals 
— W/2 there. The way 9 is modified for the Coulomb interaction reflects the change in the 
electrochemical potential of one edge due to the electrostatic potential induced by the other. 
Inserting expressions (^) in the CSLG equations (la) - (|T7J) one finds that for the short 



range interaction they are satisfied up to exponentially small terms if W ^> a -1 , where a 
is the decay constant of the single edge static solution defined previously. In the Coulomb 
case the set (||) is a static solution accurate up to terms of the order Aclxi^l/W which 
we assume to be small . The quantization of the Hall conductance is seen by integrating 
Eq. ( |Trj| ) with the result that the total current through the bar is / = a (— ¥■) ~ a o(^~) = 
Hi — — 2X c (x2 — xi) In L/W = Ji\ — fi2 + 2X c (x2 — x\) In W = v(e 2 / h)Vft a \\ where in the last 
step we restored the units of dimensions. In the Coulomb case the logarithmic dependence 
on the width of the bar indicates that the current flows predominantly in the bulk. 

We note that the sum and the difference of the values of X\ and Xi determine the charge 
of and the current through the sample respectively. In particular x\ = —X2 corresponds to 
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the physically relevant neutral Hall bar. For the set of symmetric solutions with x\ = x 2 
the edge currents balance each other while the asymmetric solutions (x\ ^ x 2 ) represent a 
system through which a net current is flowing. 

Turning to the edge excitations in the Hall bar we can now separately use the derivatives 
of the static solution (Q) with respect to either x\ or x 2 to attempt forming gapless modes 
a la Eq. @) which will be concentrated on one or another edge of the bar and propagating 
each in an opposite direction. However this prescription fails to satisfy the RPA equations 
in the case of the Coulomb interaction. Integrating the linearized continuity equation in 
this case gives the current in the x direction (0) which does not vanish on the right edge. 
We will now demonstrate that this problem is cured if the zero mode fields Sp^, Sv^ and 
5clq^ in the ansatz (f|) are taken as linear combinations d Xl p + (3d X2 p, etc., of the derivatives 
of the static solution @. This ansatz will now describe coupled density waves on opposite 
edges propagating in the same direction. Substituting such combinations into Eq. ([3d]) and 
assuming k <C W~ x we find to first order in k the dispersion relation which depends on the 
mixing coefficient /3, 



uj dp 



- + 2A, 



;i + /5)ln f— — j -(3\nW 



(9) 



k dxi 

The current pv x , Eq. (|7]) vanishes on the left wall by construction. The value of (3 is found 
by demanding that it will also vanish on the right wall. This condition leads to a quadratic 
equation in (3 giving the solutions 

and corresponding dispersion relations 



(3± = -Z±{Z 2 -l) 1 * with Z= V |fc| ' ; Ac _ V f 1 (10) 



In the case of the short range interaction, A c — > 0, these solutions tend to /3 = and (3 = — oo 
corresponding to modes concentrated on either the left or the right edge. The dispersion 
relations of these modes are the expected uj = — (dpi/ dx\)k and u = (dp 2 /dx2)k. As 

8 



anticipated the long range Coulomb forces result in inter-edge interaction which produces 
eigenmodes living simultaneously on both edges. This phenomenon has been predicted in 
Refs. [0,0 within the framework of the Luttinger model of the edge excitations in the 
QHE. It can be interpreted as a Coulomb drag of charges on one edge by the traveling 
density fluctuations on the other. Indeed @± are negative showing that the edges oscillate 
out of phase. The effect is the strongest for k — > when (3± tend to — 1 ± O (J \nk\~^. 
For a given mode the direction of propagation of the excitation is determined by the edge 
with the larger amplitude of charge fluctuations while the presence of the other edge reduces 
the velocity. For k increasing towards W~ x f3 + decreases in magnitude while (3- = f3+ l 
changes correspondingly in an opposite manner. Simultaneously the frequencies of the modes 
approach their single edge values indicating a decoupling of the edges in this limit. 

For a neutral Hall bar sustaining a small current (i.e. small x%) the current is given 
approximately by / ~ 2x\(dfli/dxi) Xl= o — 4A c £i In W. In this region we find numerically 
that dfxx/dxi ~ a + bxi where a ~ — 1.5 + 0.75(A C + A 2 ) and b is a constant close to unity over 
the range of A c presented in Fig. ^|. These observations enable us to rewrite the dispersion 
relations (|TTJ) in the following way 

bin;- 1 I 



u± = ujJ k < 



±2A 



2X c \n(W/l) -a eu c c \[ \ \k\l J 2A C 



In 



2e"T\ a 



2 



In 2 



2e"^ 
Jk\W / 



(12) 



showing that when there is no current flowing in the bar the two modes travel in opposite 
directions with equal velocities. For a nonzero total current this symmetry is broken and 
one of the modes is carried with the flow while the other is retarded by it. The dependence 
of the dispersion relation on the current suggests the possibility of a parametric excitation of 
these edge waves by driving an alternating current through the sample. For a FQHE sample 
with a length of 1mm the typical frequency of these modes is of the order of a few GHz. 
However by applying an additional constant current one can expect to excite the retarded 
branch at much lower frequencies. For the type of samples discussed here the values of the 
current for which the two terms in (|i~2"D approximately cancel each other lie in the /iA region. 
We would like to thank A. M. Finkel'stein and Y. Oreg for many useful discussions. 

9 



REFERENCES 

[1] B. I. Halperin, Phys. Rev. B 25, 2185 (1982). 

[2] M. Biittiker, Phys. Rev. B 38, 9375 (1988). 

[3] C. W. J. Beenakker, Phys. Rev. Lett. 64, 216 (1990). 

[4] D. B. Chklovskii, B. I. Shklovskii, and L. I. Glazman, Phys. Rev. B 46, 4026 (1992). 

[5] X. G. Wen, Int. J. Mod. Phys. B 6, 1711 (1992). 

[6] C. de C. Chamon and X. G. Wen, Phys. Rev. Lett. 70, 2605 (1993). 

[7] C. L. Kane and Matthew P. A. Fisher, Phys. Rev. Lett. 72, 724 (1994). 

[8] S. C. Zhang, H. Hansson, and S. Kivelson, Phys. Rev. Lett. 62, 82 (1989). 

[9] N. Nagaosa and M. Kohmoto, in Correlation Effects in Low-Dimensional Systems, eds. 
A. Okiji and N. Kawakami ( Springer- Verlag, 1994) pl68. 

[10] Yuval Oreg and Alexander M. Finkel'stein, Phys. Rev. Lett. 74, 3668 (1995). 

[11] Z. F. Ezawa, and A. Iwazaki, J. Phys. Soc. Jap. 61, 4133 (1992). 

[12] S. Giovanazzi, L. Pitaevskii, and S. Stringari, Phys. Rev. Lett. 72, 3230 (1994). 

[13] V. A. Volkov and S. A. Mikhailov, JETP Lett. 42, 556 (1985). 

[14] X. G. Wen, Phis. Rev. B 44, 5708 (1991). 



10 



FIGURES 

FIG. 1. Density (solid lines) and current density (dashed lines) profiles for the case of short 
range interactions with A s = 1. The highest curves correspond to xo = —2. Consecutive curves 
differ by Axo = 1. Note the reversal in the direction of the current when xo is increased. This 
behavior manifests a transition from a current originating from skipping orbits to a residual current 
due to uncompensated circular orbits on the rim of the Hall drop. 

FIG. 2. The chemical potential at the edge as a function of xo for short range interactions with 
A s = 1 (solid line). For the Coulomb interaction we present jl = /x + 2A c xo hiL in the realistic case 
A c = 0.03 (dashed line) and for comparison also the cases A c = 1 (dashed-dotted line) and A c = 3 
(dotted line). 
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